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Abstract 



Global perturbative QCD analyses, based on large data sets from e-p and hadron collider experiments, provide 
tight constraints on the parton distribution function (PDF) in the proton. The extension of these analyses to nuclear 
parton distributions (nPDF) has attracted much interest in recent years. nPDFs are needed as benchmarks for the char- 
acterization of hot QCD matter in nucleus-nucleus collisions, and attract further interest since they may show novel 
signatures of non-linear density-dependent QCD evolution. However, it is not known from first principles whether 
the factorization of long-range phenomena into process-independent parton distribution, which underlies global PDF 
extractions for the proton, extends to nuclear effects. As a consequence, assessing the reliability of nPDFs for bench- 
mark calculations goes beyond testing the numerical accuracy of their extraction and requires phenomenological tests 
of the factorization assumption. Here we argue that a proton-nucleus collision programme at the LHC, including a 
rapidity scan, would provide a set of measurements allowing for unprecedented tests of the factorization assumption 
underlying global nPDF fits. 
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1. The collinear factorization approach to single inclusive particle production and the role of parton distribu- 
tion functions 

The factorized QCD ansatz to proton-proton collisions establishes that the cross section for the production of a 
hadron is such collisions is given by the convolution of the incoming parton distribution functions with the hard cross 
section and with the fragmentation function. While it has been confirmed in proton-proton, when it comes to hadron- 
nucleus collisions collinear factorization is a working assumption yet to be tested. Within such approach, the nuclear 
dependence of the cross section for hadro-production, typically characterized by the nuclear modification factor 




(1) 



resides solely in the nuclear modified parton distribution functions (nPDFs) by construction (see eq.(l) in fP]). 
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Parton distribution functions play a central role in the study of high energy collisions involving hadronic projectiles 
and from the collinear factorization perspective their nuclear modifications are entirely responsible for the deviations 
from unity of the nuclear modification factor ([TJ. For protons, sets of coUinearly factorized universal PDFs are 
obtained in global pQCD analyses based on data from DIS and DY production, as well as W/Z and jet production at 
hadron colliders. On the other hand, the understanding of parton distribution functions in nuclei of nucleon number 
A, fi/A(x, Q^), is much less mature. 

It is customary to characterize nuclear effects in the parton distribution functions by the ratios 



which show characteristic deviations from unity in global nPDF analysis for all scales of tested so far. These 
effects are typically referred to as nuclear shadowing (x < 0.01), anti-shadowing (0.01 < x < 0.2), EMC effect 
(0.2 < X < 0.7) and Fermi motion (x > 0.7). 

Paralleling the determination of proton PDFs, several global QCD analyses of nPDFs have been made within the 
last decade 10, B 0. H lid based, up until recently, solely on fixed-target nuclear DIS and DY data which, compared 
to the data constraining proton PDFs, are of lower precision and cover a much more limited range of and x. Such 
limitation manifests in nPDFs with huge uncertainties and in a particularly poor gluon distribution function, since this 
can only be obtained from logarithmic -evolution, for which a wide -range is mandatory and not available so far. 

To improve on the insufficiency of nuclear DIS and DY data, recent global nPDF analyses j2|] have included for 
the first time data from inclusive high-/?r hadron production in hadron-nucleus scattering measured at RHIC iHIEIsl] 
under the collinear factorization assumption. 

In view of the importance of nPDFs for characterizing benchmark processes in heavy ion collisions, it is thus 
desirable to look for both stringent phenomenological tests of the working assumption of global nPDF fits and ex- 
perimental observables that can be used to constrain nPDFs before the construction of DIS experimental facilities 
as the LHeC [10] or the EIC ill ill . To this end we calculate the single inclusive tt" production at both RHIC and the 
LHC within the collinear factorization approach, and we ar gue that a program of hadron-nucleus collisions at the LHC 
would provide for such tests with unprecedented quality lfll ll2ll . All the calculations use LO PDFs from CTEQ6L ifisll 
with nuclear modifications EPS09LO |2] and the KKP fragmentation functions 1 14]. When the collinear factorization 
ansatz is adopted, the pj- -dependence of R'^^^ traces the jc-dependence of nPDFs. Although the precise kinematics is 
complicated by the convolution of the distributions, the qualitative dependence of the Bjorken-x of partons inside the 
nucleus with pr of the final hadron and rapidity reads x ~ pje^^ I V^- 

2. Testing collinear factorization with pA collisions 

Fig. [T] shows the nuclear modification factor for the production of neutral pions in ^[s^^ - 200 GeV 
deuteron-gold collisions at RHIC, calculated within the factorized ansatz at leading order (LO) (also consistent with 
the NLO-calculation in The RHIC data |7] in Fig. [T]have been used in constraining the nPDF analysis EPS09 
12] but they were not employed in a closely related nPDF fit |4], which provides an equally satisfactory description 
of these RHIC data. Therefore, the agreement of data and calculation in Fig. [T]is in support of collinear factoriza- 
tion: the enhancement of R^^^^ in the region around pr - 4- GeV at mid-rapidity tests momentum fractions in the 
anti-shadowing region. 

However, qualitatively different explanations of the /?J5^„ measured at RHIC are conceivable. While the above 
calculation accounted for R"i^Jjl - 0) in terms of a nuclear modification of the longitudinal parton momentum 
distribution only, such nuclear effects can be accommodated within a variety of assumptions, namely initial state 
parton multiple scattering (kr broadening) [ 15] and non-linear low-jc QCD evolution [16]. 

Yet, RHIC data being inconclusive about the dynamical explanation of the nuclear modification of the spectrum, 
we argue that repeating the same measurement at the LHC at a much higher energy and in the much wider available 



range, will help to disentangle which of the theories is correct, or in which kinematical range each one is valid 

When the single inclusive tt" spectrum is calculated at the TeV scale, the shape of the nuclear modification factor 
that would be measured at LHC mid-rapidity, Fig.|2l turns out to be qualitatively different from that observed at RHIC: 
at more than 40 times higher center of mass energy, final state hadrons at the same transverse momentum test (9(40) 
times smaller momentum fractions x,. Even tough a shift of the maximum of R" piJ^r] - 0) to such high values of pj is 




(2) 
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Figure 1: Nuclear modification factor (T) at RHIC mid-rapidity. 
Data points are PHENIX 17]. Solid line represents the EPS09 LO 
calculation. The dashed lines are the EPS09 uncertainty (inner) 
plus data normalization uncertainty (outer). 



Figure 2: EPS09 LO prediction for the nuclear modification factor 
of neutral pion single inclusive production in proton-lead collisions 
at LHC mid-rapidity. All LHC rapidities are in the center-of-mass 
(shift to lab Av ~ 0.47). 



a natural consequence in collinear factorization reflecting nuclear modifications of PDFs, it cannot be accounted for 
in terms of A:7-broadening, since such approaches predict a mild dependence of the spectrum with the center-of-mass 
energy |15]. 

Moreover, when the single-inclusive spectrum is calculated for different values of rapidity the results. Fig. [3] show 
that the rapidity dependence of R^pj^ allows one to scan the main qualitatively different ranges of standard nPDFs in 
an unprecedented way. The behavior of the nuclear modification factor when saturation is taken into account f|l6i1 is 
rather different: a strong suppression factor is found which is incompatible with any of the existing nuclear parton 
distribution sets. 
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Figure 3: Rapidity dependence of S^pj Q for y.VNN = 8.8 TeV pPb (LHC). The different plots scan the dependence from v = -4 (close 
to Pb projectile rapidity) up to y = 4 (close to proton projectile rapidity). Labels indicate whether the nuclear modification originates mainly 
from the shadowing (S), anti-shadowing (AS) or EMC regime. Vertical lines illustrate the rapidity-dependent pr range, which can be accessed 
experimentally with more than Wevcnts = 1000 (= 10) per GeV-bin within one month of running at nominal luminosity. 

Thus, a measurement at the LHC which shows the rapidity dependent behavior in Fig. |2] and [3] would be in strong 
support of collinear factorization as opposed to initial state multiple scattering and CGC approaches. 



3. The LHC as a tool to constrain the gluon nPDF 

Despite the perspectives for qualitative tests of collinear factorization, current global analyses of nuclear parton 
distribution functions come with significant uncertainties as we can see in Fig. 12 in We have compared in Fig.|4] 
the nuclear modification factor for two nPDF sets, which show marked differences: in contrast to EPS09 the gluon 
distribution of HKN07 |0] does not show an anti-shadowing peak but turns for x > 0.2 from suppression to strong 
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enhancement at initial scale = I GeV^. Fig. |4] thus illustrates that within the validity of a coUinearly factorized 
approach, LHC data can resolve the qualitative differences between existing nPDF analyses and can improve signifi- 
cantly and within a nominally perturbative regime on our knowledge of nuclear gluon distribution functions. 
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